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Abstract 

We test the AdS/CFT correspondence by computing the partition function of some N = 2 quiver 
Chern-Simons-matter theories on three-sphere. The M-theory backgrounds are of the Freund- Rubin 
type with the seven-dimensional internal space given as Sasaki-Einstein manifolds Q 1 ' 1,1 or V 5,2 . 
Localization technique reduces the exact path integral to a matrix model, and we study the large- iV 
behavior of the partition function. For simplicity we consider only non-chiral models which have 
a real-valued partition function. The result is in full agreement with the prediction of the gravity 
duals, i.e. the free energy is proportional to iV 3//2 and the coefficient matches correctly the volume 
of Q 1 ' 1 - 1 and V 5 > 2 . 



PACS numbers: 11.15-q,11.25.-w 



I. INTRODUCTION 



We have seen remarkable progress made about M2-brane dynamics recently. It is now 
well established that M = 6 Chern-Simons gauge field theory (ABJM model) with U(N)k x 
U(N)-k correctly describes M2-branes at orbifold singularity C 4 /Zfc [Ij. But as usual with 
AdS/CFT correspondence, many of the checks with maximally possible duality relation 
strongly depend on the symmetry, and for a more convincing and challenging test, one 
would like to study models with less supersymmetries. 

Indeed, there have been a number of works which put forward new, more nontrivial 
duality relations regarding M2-branes [2"HT2"]. The richest and still challenging category is 
probably Af = 2 models. The dual geometry is usually of Freund-Rubin type, i.e. AdS4 x X 7 
where X 7 is a seven-dimensional manifold satisfying the Sasaki-Einstein property. One can 
use the machinary of toric geometry [13], brane tiling [13] , crystals [15] etc., but we are yet 
to have a complete understanding even with the more tractable class of toric Sasaki-Einstein 
manifolds. Most of the earlier works just checked that (a particular branch of) the vacuum 
moduli space coincides with the gravity side, and studied the chiral ring structure perhaps 
including some nonperturbative states involving baryonic operators. 

The situation has changed recently, mainly due to important progresses on the field 
theory side. One is the computation of supersymmetric indices, where one can certainly 
check beyond the chiral ring. Indeed, some of the models which showed precise vacuum 
moduli space turned out to give wrong indices and are ruled out. For more details, see |16l - 
[22] . The other route is the computation of exact partition function on three-sphere, and it 
is the aim of this paper to present calculation of the partition function for some AdS4/CFT3 
proposals with four supercharges. 

It has been known for a long time that in M2-brane dynamics the number of degrees of 
freedom scale as A 3//2 , where N is the number of M2-branes. First of all, using localization 
technique it is shown that the partition function on three-sphere can be expressed as a 
matrix integral [23J. Using an exact resolvent it was first confirmed in a beautiful paper 
[21] that the free energy has the right behavior of A 3//2 . For more details and related works, 
readers are referred to [251428] . 

For more general models, we do not have the luxury of an exact resolvent, but a powerful 
numerical and analytic method have been proposed in [2H]- In that procedure, the matrix 
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integral is approximated using saddle point method. The saddle point equations are turned 
into an auxiliary first-order coupled differential equations which quickly give the information 
on the behavior of roots. In [2S] the authors were interested in tri-Sasakian geometry and 
J\f > 3 quiver Chern-Simons models. For the models studied in j29] , the free energy exhibits 
nice N 3 / 2 behavior and the coefficient is consistent with the volume of seven-dimensional 
internal space through the well-known AdS / CFT dictionary 



In these models with J\f > 3, all the fields are equipped with canonical conformal di- 
mension and R-charge, i.e. A — R — 1/2. But in A^ = 2 theories in general the fields 
can acquire anomalous dimensions. The problem of putting general superconformal M = 2 
theories on three-sphere and computing the partition function has been solved in [301 E] ■ 
More concretely, it is shown that chiral multiplets with an arbitrary R-charge can couple to 
gauge fields and curvature of S 3 in a way preserving OSp(2\2) x SU(2). This symmetry is 
a half of the superconformal group in addition to the isometry of S 3 , SU(2) x SU(2). 

Quoting eq.(1.2) of [30J, the partition function is given as 

Z = J \\ du e^ Tru2 det Ad (sinh(7™)) J] det fli (^ (1 - A * +m) ), (2) 

Cartan Ri 

where Aj is the conformal dimension of i-th chiral multiplet in gauge group representation 
Ri. The function £(z) is known as double sine function, see for more details (32j 133] . 

One can in principle try to compute Z for all the Chern-Simons dual proposals. When 
we started working with them, we have found that chiral-models which give complex-valued 
Z are rather difficult to deal with. We do not have the full understanding of how the roots 
condense in large-iV limit yet. On the other hand, it turns out that simpler non-chiral 
models with two-nodes produce dependable but still nontrivial results. In other words, we 
can verify the free energy agrees with ([I]) in large- N limit. 

In this paper we study three AdS4/CFT 3 examples. In Secjlljwe study the first example 



which is dual to V 5,2 , where the gauge theory comes with adjoint representation. In Sec HI 
we study models involving chiral multiplets in fundamental representations. One is dual 
to Q 1 ' 1 ' 1 , and the last example provides an alternative description of V 5,2 . We give a brief 
discussion on our results in SecJfVl 
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FIG. 1: The quiver diagram for Chern-Simons dual of AdSi x V 5,2 

Note added: There are two papers on the arxiv which have significant overlap with this 
article. As this paper was being typed, we became aware of a work [31] and we coordinated 
the release of the paper. Later we received [35], which provides a comprehensive analysis 
for non-chiral Af = 2 models. In particular, the F-theorem in [35] clarifies a misconception 
on Z-extremization in the original version of this paper. 

II. QUIVER WITH ADJOINT MATTERS: V 5 ' 2 

In this section we consider the dual field theory for V 5 ' 2 geometry. This geometry is a 
direct higher-dimensional generalization of conifold. The cone over V 5,2 is a singular space 
defined as 

z\ + z 2 + z% + zt + z 2 = 0, Zi £ C. (3) 

Obviously this geometry has U(1)r x £0(5) isometry. V 5,2 is an example of non-toric Sasaki- 
Einstein manifolds. According to AdS/CFT, M-theory in AdS4 x V" 5 ' 2 is dual to a strongly 
coupled three-dimensional M = 2 conformal field theory. 

The dual as a quiver Chern-Simons system is proposed in [10J, see FigJT] The gauge 
group is U(N) x U(N) with Chern-Simons levels (k,—k). The matter fields include six 
chiral multiplets in total: A a , a = 1, 2 in (N, N) representation, B a , a = 1, 2 in (jV, N), and 
we also have one chiral multiplet in the adjoint representation for each gauge group, called 
$i,$2- The vacuum moduli space provides an abelian orbifold V 5,2 when one chooses the 
following cubic superpotential 

W = Tr($? + $3 + < i 'i<\> : .\,I>> l + QzBiAj)). (4) 

For a study of this duality relation using rotating membrane solutions, see Ref . [36J . 

It is straightforward to write down the formula for the partition function using the recipe 
given in [3TJI I3"T] . with arbitrary R-charges or equivalently the conformal dimensions of the 
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matter fields. From the symmetry of our system, it is obvious that the bifundamental fields 
have the same R-charges, and we will call it A&. The two adjoint fields $1, $2 should also 
come with the same R-charge A a . Now if we designate a Cartan subgroup for U(N) x 
U(N) using eigenvalues Aj, Xi, i = 1, ■ ■ ■ ,N, the partition function is given by the following 
integral. 

2 
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The function £(z) was used in [30| l3Tj to express a one-loop determinant for a chiral 
multiplet with a generic value of R-charge. As an infinite product, it is defined as 



exp t[z ) 



More concretely it can be expressed as [30 



n 

71=1 



n + z 



n — z 



(6) 



^) = ln(l - e 2 " 2 ) + - [nz 2 + -Li 2 (e 

2 V 7T 



Z7T 
12' 



(7) 



Although it is not immediately seen from this expression, from the definition as infinite 
product it is clear £ is an odd function, and also £(z) 6 K if z G R. 

For the manipulation of the integrand for Z, we find the following relation of £(z) very 
useful. If z = x + iy with y £ K, we can derive 



(x + iy) + £{x — iy) = £{x + iy) — £{—x + iy) 



-2nx\y\ + e 



■27rn|s/| 



n=l 



2a; cos(27tox) 
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2|y| 1 
— + — 



7? 



7T7?/ 



sin(27rnx) 



(8) 



Let us check the convergence of multiple indefinite integrals in ([5]). First of all, the Chern- 
Simons part with k dependence is highly oscillatory, but in itself the integral is convergent, 
since for instance sin x 2 = a/vt/2. Recall that we are doing the integration along real- 
axis. The remaining part in the integrand with sinh is exponentially growing, and the 
part with £ is exponentially decaying, as we send A or A to infinity. Suppose we take one 
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particular eigenvalue A» to infinity, and keep others finite. Then the integrand, apart from 
the oscillating part, behaves like 



exp[(iV-l)(A a + 2A b -2)A i ] (9) 

Then it is obvious that the integration in ^ is divergent if A a +2A^ > 2. When A a +2A& < 2 
it is absolutely convergent, and we expect Z becomes smaller as we make R-charges smaller 
and smaller. We are interested in the marginal case of A a + 2A& = 2 in this paper, but will 
make some comments about A a + 2A& < 2 case later. 

We are interested in the behavior of Z in the leading order as iV — > oo. We will employ 
the saddle point approximation and take the large- N limit, closely following the procedure 
described in [22] • The saddle points are determined by setting the following quantities to 
zero. In the following F is related to Z by Z = J e~ F . 



dF ik . v-^ Xi - Xj 

— — = Xi — > coth 

dXi 2vr ^ 



+ H 1 - A *> ^r 1 ) + m ^ - A - ^r 1 )] > ( 10 ) 



j 

\ . _ \ . 

•j 



dF ik ~ ^-^ Xi — X 
Xi — c °th 



8Xi 2tt ^ 2 



- [2m(l - A 6 , ^— ^) - m(l - A a , ^— . (11) 



A V^) " -(1 - A a , 

3 

In order to avoid too much clutter we introduced a shorthand notation 

i 

m(x,y) = -— (£'(x + iy) -£'(x-iy)) 
y sin 2irx — x sinh 2ny 



(12) 



cos 2nx — cosh 2ny 
where we used 

£'(z) = -TTZCOt(lTz). (13) 

From the saddle point equations above, one easily notices that the saddle point equations 
are invariant under (i) (Aj,Aj) — > — (Aj,Aj), and (ii) (Aj,Aj) — > (A*, A*). The implication is 
that the root distribution is symmetric with respect to the origin, and Xi is the complex 
conjugate of Aj. 
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N=20 and 100, k=l, A„=Aj= 2/3 
Im(A) 



N=20 and 100, k=l, A„=l/2, A*= 3/4 
Im(A) 
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Re(A) 



FIG. 2: These two figures show saddle points for V 5,2 theory with different R charges, from our 
numerical analysis. The roots are plotted with black and blue dots for N = 20 and 100, respectively. 



which are represented by green lines. 

A practical method to find the saddle points is to treat the eigenvalues Aj, Aj as functions 
of time t, and solve an auxiliary coupled ordinary differential equations (ODE) 



By choosing the initial conditions Aj(0), Aj(0) appropriately, the eigenvalues as functions of 
t converge rapidly as t — > oo. Obviously the fixed points Aj(oo), Aj(oo) would give us the 
roots of saddle point equation. This method was utilized by the authors of Ref . to study 
the large- A" behavior of partition function for M = 3 Chern-Simons quiver theories. We find 
this technique is dependable for the non-chiral models we deal with in this paper, and we 
closely follow the procedure in [29J. 




We put the coupled ODE on computer and have obtained the roots. They condense and 
make a cut, and the roots for some different values of R-charges are displayed on complex 
plane in Figjij We first note that the distribution is symmetric under (Aj, Aj) — > — (Aj, Aj), 
and in fact it always holds that A* = Aj. This is of course not unexpected, considering the 
symmetry of the integral mentioned earlier. 

As we analyze the numerical data for different values of N, we find that the real part of 
A scales as N a with a ~ 1/2, and the imaginary part stays within the same interval of order 
1. This behavior is qualitatively the same as the case of ABJM model reported in |29j . and 



The real value of the roots scales as \flSf . The imaginary of part of the roots is bounded by ±Afr7r 



d\i _ dF d\i _ dF 



(14) 




7 



we also write 



\ i = N a Xi + iyi 



(15) 



For convenience, we assume we have re-arranged the roots so that Xi < Xj if i < j. As N 
goes to infinity, Xi, yi become continuously distributed over an interval. 

Another notable feature of Fig|2]is that the range of Im(Aj) changes for different values 



of R-charge. This can be understood if we go back to the saddle point equations (10 ),(!!) 



Let us look at the following term from the second line of (10). 

Xi — Xj 



E 

j 



m(l - A 6 , 



2tt 
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With (15), for j — % this expression contains 

% y sin(27r(l - At)) - tt(1 - A b ) sin(2y) 

m l - A 6 , i Vi it) = — — — — , 

7r cos(27r(l — At)) — cos(2y) 

which is divergent at y — 7rAt. At this point the roots feel infinite force, and the roots 
cannot cross this boundary. This argument is verified in Fig|2} 

Now we take the large- iV limit of F analytically. We first write F = F ext + F int , where 
F ext describes the Gaussian part with a single summation over the eigenvalues, and F int 
describes the interaction of the eigenvalues with a double summation. Juggling with the 
summation indices, we can rewrite Fi nt in the following way. 
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21n(l - e 
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2vr 7 
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+2£(1 - At + z 



2tt ' 
2vr ' 



2£(-l + At + 2 
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2tt ' 
Aj — Xi , 



2tt 



+ - A.) + 2£(1 - At + *^^) + 2*(1 - At - i^F 1 )}- (18) 

i 

It turns out that the terms within the double summation give 0(N 2 ~ a ) contribution, 



and the single summation part in the last line of (18) is sub leading as 0{N). Now for the 



terms involving the function t for lines 2-5 of (18), we can make use of the relation rt8 



One discovers that the linear terms such as In e Aj Al + In e Xj Xi in the first line get exactly 
canceled by similar terms from the matter fields, if we impose 



A a + 2A b = 2. 



(19) 



This is essentially the same asymptotic property we discussed around Eq.([9]). It is certainly 
consistent with the superpotential Q, since Ty(ABQ) should have R-charge 2. In fact, 
considering everything in (J3J) one easily sees we should set A a = A b = 2/3, if we are to 



relate this system to V 5 ' 2 . But we will just assume (19) and leave A b undetermined in the 
following computation. 

Now the rest of the summand is given as a series expansion, and we move on to take 
large- N limit. We introduce a continuous variable s instead of summation indices i,j, 



s = i/N, < s < 1. 



(20) 



Then a summation is transformed into an integral, 



where = Max(xj) and p(x) 



1 - r 1 

i=i ^° 



F(s(x))p(x)dx, 



(21) 



ds 
dx ' 



Ignoring the terms in the last line of (18), we obtain 
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(22) 



where we defined a shorthand notation a = 27r(l — A a ), (p b = 27r(l — A b ). 



Eq.(22) looks complicated, but if we integrate by parts repeatedly we can obtain an 
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expansion in terms of 1/N a . Consider for instance 
dx'e- n ^- x ^F(x') 

fai e -nN a (x-x') e -in(y(x)-y(x')) J?( x ') 
_ X \ c -nN a (x-x')-in(y(x)-y(x')) p( x ')] +0(1/N 2a ) 

nN a L V J \x*=-x* V ' ; 

= -J^F(x) + 0(l/N 2 «). (23) 

The above approximation is valid if in F(x), x's do not appear with a factor of N a . For 
instance the density function p(x) satisfies this criterion. In this way we can remove the 
integration over x' in (22), and arrive at an expression 

F int = N 2 ~ a dxp(x) 2 f(2y(x)), (24) 

J — a;* 

where the computation leads to 

1 



oo r 

f(2y(x))=Aj2\-, 

n=l I 



cos(n0 a ) H r sin(n0 a J 



27in 2 Tin 3 

— cos(n0 b ) cos(2ny) H ^— sin(n0 b ) sin(2ny) H sin(n0b) cos(2ny) (25) 

Tin 1 Tin 1 -nn 6 

One can sum this expression using Fourier series of simple linear or quadratic functions. 
The answer can have a discontinuity at y = <pb/2 in principle, but rather miraculously, the 
discontinuities cancel out and we obtain the following simple result. 

f(2y) = -0 fe ((27T - b ) 2 - % 2 ), -7T < y < tt. (26) 

7T 

Apart from the fact that the intercept at y = is a function of R-charge assignment, f(2y) 
being a quadratic function is exactly the same as the result of ABJM model in large- limit 
reported in [29|. 

For the extremization of the whole action we also consider the external force described 
by the Gaussian part. It is given as 

h N 

8=1 

kN 1+a f x * 

/ dxxp(x)y(x) (27) 
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FIG. 3: These figures show the density of roots p(x) for V 5 ' 2 . The dots are numerical data and 
the red lines show the analytical predictions. 

Comparing F ext against Fi nt , it is clear the sum should be extremized when a = 1/2, which 
is consistent with numerical analysis. 

Taking into account the constraint that the total number of roots is N, we have the 
following expression 



jV 3 /2 



k 



dxxp(x)y(x) + / dxp(x) f(2y(x)) 

IX I / 27T 



dxp(x) — 1 



(28) 



where \l is a Lagrange multiplier. In fact this is identical to Eq.(16) of [29], and holds in 
general for any non-chiral models with two nodes we will consider in this paper. It is easy 
to solve the equation of motion, and one finds the density function is constant while y(x) is 
linear, i.e. 



p(x) 



, A;(2vr-0 b ) 2 
y{x) = x 



(29) 



4</» 6 (27T - <j) h f ' aK ' 2/2 

One can check if p(x) is really a piece-wise constant function or not, with numerical data 
It is indeed the case, as one can see in Fig|3j 
Then F is expressed as a function of x*, 



F = N 3/2 



<P b (2* ~ <p b ?_ + k 2 xl 



27TX* ' 247T06 

The extremization point is easily found, and with the corresponding imaginary part it is 

1 
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.1 jJs 
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k ' °* 2' 

Here we see that Max(Im(y)) is indeed bounded by 7rA{,, in consistent with our wall of 



infinite force argument around Eq.(17) 
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Finally the answer for free energy is 

F = w/a #Fjj! (32) 

3tt V ; 
Now if we substitute the natural value A a — A^ — 2/3, the free energy becomes 

' 16 V/ 2 iV 3 / 2 , (33) 



27 

and it is certainly consistent with the general formula with the volume of V 5,2 given in [10 



Vol(v5,2) = w- (34) 

III. QUIVER WITH FUNDAMENTAL MATTER: Q 1 1 1 AND V 5 ' 2 

A. Q 1 ' 1 ' 1 

Our second example is a toric homogeneous Sasaki-Einstein manifold Q 1,1,1 . This space 
is a £7(1) fibration over CP 1 x CP 1 x CP 1 , so the isometry group is U(1) R x SU(2) 3 . There 
are several proposals for the field theory duals of AdS4 x Q ' ' background. A four-node 
quiver is proposed for a specific abelian orbifold Q 1 ' 1 ' 1 /Zk in Ref. [37], which also studied 
a closely related geometry, Q 2 ' 2,2 = Q 1,1 ' 1 /^. These models are chiral, and although it 
is straightforward to write the integral formula for their partition functions, we found it 
difficult to perform both numerical and analytic computations. 

We consider instead a non-chiral model proposed in [TTl [T2"] . Like our previous example, 
this theory is closely related to ABJM model. In addition to the bifundamental fields A4, B^ 
one adds a pair of ( ant i) -fundamental representation fields to each gauge group node. The 
quiver diagram is given in Fig|4| In addition to the ordinary quartic superpotential of ABJM 
model, one adds cubic terms 

W = e ij e kl TrAiB k AjB[ + Q 1 A 1 Q 1 + Q 2 A 2 Q 2 . (35) 

It is argued that this theory gives C(Q 1 ' 1 ' 1 ) when the bare Chern-Simons levels vanish. 
Unlike our previous example V" 5 ' 2 in the last section, even if we are using the fact R(W) = 2 
we cannot determine the R-charges completely. Using the SU(2) x SU(2) symmetry, we 
may set R(Q { ) = A f ,R(Qi) = Af, D(A { ) = A A , and D(B { ) = A B . Using the fact W 
is a marginal operator, we may expect Af + Aj + A a = 2A^ + 2A B = 2. If we assume 
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FIG. 4: The quiver diagram for Chern-Simons dual of AdS$ x Q 



1.1,1 



the bifundamental fields retain the canonical conformal dimension 1/2, we would get A A 
A B = 1/2, A/ = 3/4. 

We can write down the partition function as follows. 
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x J^J exp 
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Z7T Z7T 



2€(1- A/-i-i) + 2£(l- A f - + 
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f 2,71 ' 



One can easily check this expression is real- valued, and satisfy the following property. 



(36) 



Z{A A ,A B ,A f ,A f ) = Z{A A ,A B ,A f ,A f ). 



(37) 



Let us now check the convergence of this integral. As we take one of Aj to infinity while 
keeping all others finite, we see the leading behavior is 



exp [{(N - 1)(A A + A B - 1) - 2(1 - A f )} A, 



(38) 



Then clearly for large N the integral is divergent for A A + A B > 1 . A A + A B = 1 is the 
marginal case we are interested in, and then the convergence depends on the value of Af: 
We have convergence if Af < 1. And it is the same with A and Af. 

With the saddle point equation, one can check that A* = A» if we assume Af = Af. In 
order to simplify the following analysis, we will assume this is the case from now on. The 
root distribution from the numerical results are given in Fig|5j 
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FIG. 5: The distribution of roots for Q ' ' theory with different R charges. The roots are plotted 
with black and blue dots for N = 20 and 100, respectively. The real part of the roots scales as 
\/7V. The maximum values of imaginary part of A are bounded by itAb, which are represented by 
green lines. 
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FIG. 6: These figures show the density function p(x) for Q 1 ' 1 ' 1 . The dots are numerical data and 
the red lines show the analytical predictions. 

Now we turn to the analytic computation. Taking the continuum limit, we first note that 
the linear terms cancel if Aa + A# = 1. We assume this is true, and after some manipulation 
similar to SecJTTl we arrive at 



F = iV 3 / 2 



dx I 1 - A f - 



dxp(x) 2 f(2y(x)) + 2 
The interaction part is described by 

f(2y) = 7T 2 - (2y + 2irq) 2 
where we set for convenience A^ = 1/2 — q, A B = 1/2 + q. 



y[x) 
2tt 



p(x)\x\ 



(39) 



(40) 
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Now solving the functional variation equation, we obtain 

7r 2 \x\ a - (2 - 2A f + q)2n\x\ . . 

y= (2-2A f + q )4n\x\-2»- q *> P = 4^ (41) 

Unlike the previous example, the density function is not constant, and this analytic calcu- 
lation is supported by numerical results. See Fig|6j 

Then we fix /i from the normalization condition J p(x)dx = 1, and find the extremal 
point of F as a function of x*. Finally, the free energy as a function of Af,q is 

p = N3/2 V2n(5 -4A / + 2 9 ) 
3 v /3-2A / + g 



If we substitute A a + 2 A/ = 2, which is consistent with superpotential Eq.(35) 



F = N ^±L. (43) 
3v 3 

This matches nicely with the volume of Q 1 ' 1,1 , which can be found for instance in [38J. 

VoliQ 1 *' 1 ) = (44) 

Note that, unlike the previous example of V$ } 2, the right value of F does not uniquely 
determine the R-charge values for all matter fields. 



B. V 5 ' 2 

We already studied a Chern-Simons-matter theory which describes V 5 ' 2 in Secjll} It is 
possible to devise another dual field theory using matter fields in fundamental representation 
[TT| IT2] . One starts with pure super- Yang-Mills theory without Chern-Simons terms, with 
gauge group U(N). There are three chiral multiplets X, Y, Z in adjoint representation, and 
we consider adding k fundamentals Q a and antifundamentals Q a , with a — 1, - • • , k. 

The superpotential is given as 

k 

W = TrX [Y, Z] + Qa{X 2 + Y 2 + Z 2 )Q a . (45) 
a=i 

If we are to use this superpotential, the R-charge is obviously 2/3 for adjoint fields, and 1/3 
for (anti)-fundamental fields. The vacuum moduli space correctly reproduces C(Y 5 ' 2 )/Zfc 

mm- 
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X,Y,Z 



FIG. 7: The quiver diagram for Chern-Simons dual of AdS^ x V 5,2 



The matrix integral is written as 



1 

M 



n^n 

i=l i<j 



2 sinh 



A,; — A, 



n 



exp 



M(l-A a + i X% Xj 



2tt 



x 



n 



exp 



k£{l 



A f 



2tt' 



(46) 

R(Q). This 



where A a is the R-charge of adjoint fields X,Y,Z, and Af = R(Q),Aj> 
expression is always real-valued. 

Let us again consider the convergence of the integral. As we send one particular A, to 
infinity while keeping other variables finite, the leading contribution comes from the first 
line, and the integrand behaves 



exp [{(3A a -2){N-l)-k(2-A f - A f )/2} Aj 



(47) 



So it is obvious that for A a > 2/3, the integral in (|HJ is divergent. As A a becomes smaller, 
Z should decay faster and faster. We are basically interested in the large-iV limit with 
A a = 2/3, and perform the computation with an arbitrary R-charge for quarks, i.e. Af = 
A j = 1 — 2h. Note that, once we set A a = 2/3, the integral is convergent for h > 0. 

One can proceed in the same way as before, and from the saddle point equation it is 
clear that the roots are actually real- valued. The linear terms like lne Ai_Aj get canceled by 
similar terms from adjoint fields, since A a = 2/3. 

As the eigenvalues are real- valued at the saddle point, the free energy only contains the 
density function p, and the equation of motion in the continuum limit is much simpler than 
previous examples. From numerical computation we have Aj = N a Xi and a ~ 1/2. Analytic 
computation shows a = 1/2 in the same way as before. 

The interaction part comes from the gauge field contribution and the matter fields in 
adjoint representation. They are of double- summation form and turn into a multiple inte- 
gration, but in the large- iV limit we can remove one integration by partial integration and 
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FIG. 8: The left figure shows saddle points for flavored V 5 ' 2 theory with N = 100 and R charge 
A a = 2/3 and Af = Aj = 1/3. All roots are real. The right figure shows the density function 
p(x). The dots are numerical data and the red line shows the analytical predictions. This show 
that numerical data and analytical prediction p(x) = 3(4ir — 3|x|)/(167r 2 ) are same. 



evaluating at the boundary. The final result turns out to be 

16vr 2 iV 2 - a 



int 



27 



dxp(x) . 



(48) 



Now the quark fields Q, Q give single-summation expression in F, and in fact they play a 
similar role as the Chern-Simons terms. In the large- iV limit the leading contribution is 



fund 



2hkN 



l+a 



dxp(x)\x\. 



Now we set a = 1/2 and introduce a Lagrange multiplier to write 



dx ■•'..>- 



27 



p(x) + 2hk / dxp{x)\x 



p_ 
2tt 



From the equation of motion, the density function is obtained 

27(/i-4/iibr|x|) 



p(x) 



64vr 3 



(49) 



(50) 



(51) 



Note that it is a piece-wise linear function. This fact is confirmed by numerical results, see. 
Figj8] Now we can determine p, from the normalization of p(x). The total free energy as a 
function of is 



F = N 3/2 [ hkx* - 



9h 2 k 2 xl 8vr 2 



32vr 2 ' 27a:*. 

The final solution for free energy extremized with respect to x* is 

16tt 



9^ 



{hkf' 2 N 3 l 2 . 



(52) 



(53) 
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Now if we use the superpotential and set A/ = Aj = 1/3 or h = 1/3, this is the same as 
(33) and consistent with the volume of V 5,2 and AdS/CFT. 



IV. DISCUSSION 

We have provided a test of AdS 4 /CFT 3 correspondence for three Af = 2 models in 
this paper. The dual geometries are Q 1 ' 1,1 and V 5,2 , which are both relatively simple, 
homogeneous Sasaki-Einstein manifolds. The partition function is calculated in large-iV 
limit, and in the leading order the free energy with appropriate R-charge assignment shows 
nice agreement with the predictions of AdS/CFT. Compared with the study involving just 
the chiral ring structure, our result renders very strong support to the credibility of the dual 
field theories. 

For explicit computations, we first used symmetry argument and marginality of the matrix 
integral to reduce the number of independent R-charges, and expressed Z as a function of 
remaining .R-charges. We can then read off a condition for cancellation of linear terms from 
the integrand, like A^+Ag = 1 in the Q l ' l)1 model. If this condition does not hold, the linear 
terms with double summation survive and in fact dominate the action. For A^ + A# < 1, 
the integral is convergent and this corresponds to a rather more conventional matrix model 
integral. The roots exchange usual log \\ — \j\ type interaction, and the real parts as well 
as the imaginary parts are bounded within a finite interval, as iV — > oo. 

We have checked this prediction numerically. For the Q 1 ' 1 ' 1 model we have studied here, 
numerically fitting the data for N = 20,40, • • • , 120, we have obtained a = 0.569391 for 
A^ = A# = 1/2, and a = 0.104406 for A^ = A# = 1/4. This makes the free energy 
scale as N 2 below the line of marginality. In Figj9] we have illustrated this fact using our 
numerical data. At a point on the marg inal line, A A = A B = 1/2, F ~ N 3 / 2 as we have 
argued already. But away from it, for instance at A^ = A B = 1/4, the numerical result 
gives F ~ N 2 . If one goes over the line of marginality, with larger R-charges, then Z — > oo 
or F — > — oo as we have argued already. 

After imposing conditions on R-charges using symmetry and the marginal convergence of 
Z, we have performed analytic computation and found expressions for F/N 3 ^ 2 as a function 
of the remaining R-charges. For the examples considered in this paper, if we make full use 
of the marginality of the superpotential in the dual field theory, F/N 3 ^ 2 agrees exactly with 
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FIG. 9: These figures how free energy scales with respect to N for different R charges for Q ' > 
model. The graphs are log-log plots. The blue dots are numerical data and the red lines are fitting 
lines. We took data for N = 20, 40, 60, 80, 100, 120, 140, 160. 

the prediction of AdS/CFT Eq.Q. 

In [30], the author put forward an interesting argument that the partition function as a 
function of conformal dimensions is extremized at the correct values. In the large- N limit 
we are interested in here, Z itself converges to zero and one needs to consider a refinement of 
Z-extremization using J 7 = — limjv^ 00 (ln Z/N 3 ^ 2 ) instead. This is addressed in [35] 1 , where 
it is shown that the free energy in the large- N limit for strongly coupled M = 2 non-chiral 
quiver Chern-Simons is consistent with F-extremization theorem. 

To check F-extremization, one should consider R-charge assignments consistent with the 
marginality of superpotential, and then F should be extremized by the correct values of the 
remaining freely-adjustable R-charge values. Our result is in harmony with F-theorem. In 
fact, when we use both symmetry argument and marginality of superpotential, all R-charge 
values are either completely determined (for V 5,2 models) or F has no dependence on the 
remaining choice of R-charge values (for flavored Q 1,1 ' 1 model). 

In order to see the extremization more clearly, for instance with our first example we 
need to calculate more generally Z(Aa, A b , A $1 , A$ 2 ) just using marginality of W, i.e. 
Z(Aa, A b = 4/3 — A^, A$ 1 = 2/3, A$ 2 = 2/3). But it is obvious that this function of A^ 

1 We thank D. Jafferis who pointed out erroneous statements on Z-extremization in earlier versions of this 
paper. 
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should be extremized at Aa = 2/3, due to symmetry argument. 

In this paper we studied only non-chiral Chern-Simons quivers with real-valued partition 
function. For chiral models such as the triangular quiver for M 3 ' 2 model [3J 0], it is more 
difficult to identify the behavior of roots and perform analytic calculation if we apply the 
same method used here. We hope to be able to report definite answers for chiral model 
proposals in future publications. 
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